Tunable backaction of a dc SQUID on an integrated micromechanical resonator 
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We have measured the backaction of a dc superconducting quantum interference device (SQUID) 
position detector on an integrated 1 MHz flexural resonator. The frequency and quality factor 
of the micromechanical resonator can be tuned with bias current and applied magnetic flux. The 
backaction is caused by the Lorentz force due to the change in circulating current when the resonator 
displaces. The experimental features are reproduced by numerical calculations using the resist ively 
and capacitively shunted junction (RCSJ) model. 



It has recently been demonstrated that a macroscopic 
mechanical resonator can be put in a quantum state [l[ 
by coupling it to another quantum system. At the same 
time, linear detectors coupled to mechanical resonators 
are rapidly approaching the quantum limit on position 
detection. This limit implies that the resonator position 
cannot be measured with arbitrary accuracy, as the de- 
tector itself affects the resonator position [2]. This is 
an example of backaction. Backaction does not just im- 
pose limits, it can also work to one's advantage: Back- 
action can cool the resonator, squeeze its motion, and 
couple and synchronize multiple resonators. Different 
backaction mechanisms have been identified: When using 
optical interferometers or electronic resonant cir- 

cuits (tJ-EI, backaction results from radiation pressure. In 
single-electron transistors (SET) [lOj, Cooper-pair boxes 
carbon nanotube quantum dots [HI, [l3| , or atomic 
and quantum point contacts [13, [Hj backaction is due 
to the tunneling of electrons. Recently, we have used 
a dc SQUID as a sensitive detector of the position of 
an integrated mechanical resonator [16| . This embedded 
resonator- SQUID geometry enables the experimental re- 
alization of a growing number of theoretical proposals for 
which a good understanding of the backaction is required 
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In this Letter, we present experiments that show that 
the dc SQUID detector exerts backaction on the res- 
onator. By adjusting the bias conditions of the dc SQUID 
the frequency and damping of the mechanical resonator 
change. The backaction by the dc SQUID has a different 
origin than in the experiments mentioned above: It is due 
to the Lorentz force generated by the circulating current. 
Numerical calculations using the RCSJ model for the dc 
SQUID [24| reproduce the experimental features. 

The device (Fig. QJi) consists of a dc SQUID with 
proximity-effect-based junctions [l6j. A part of one arm 
is underetched, forming a 1 MHz flexural resonator with 
length i = 50 /im. In this Letter we present data on 
a device in an in-plane magnetic field of B = 100 mT. 
Measurements have been performed at several magnetic 
fields and on an additional device; the observed back- 
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FIG. 1: (a) Schematic overview of the dc SQUID with the 
suspended beam and measurement setup. A magnetic field B 
transduces a beam displacement u into a change in magnetic 
flux. A bias current Ib is sent through the SQUID and its 
output voltage is measured. The flux $ is fine-tuned with a 
stripline current If that is controlled by a feedback circuit 
(dashed) that keeps the output voltage V at Vsp- (b) The 
bias current dependence of the measured V m i n and Vmax- (c) 
The amplitude (bottom) and phase (top) response. The line 
is a fitted harmonic oscillator response [33j] . 



action is similar [25j. First the dc SQUID is character- 
ized. The output voltage of a dc SQUID depends on the 
magnetic flux through its loop [24] and we measure 
the minimum and maximum voltage (Vmin and F max ) by 
sweeping the flux over a few flux quanta $0 = h/2e with 
a nearby stripline (Fig. la). This is repeated for dif- 
ferent bias currents to obtain the current- voltage curves 
shown in Fig. [TJd. The maximum critical current is 
jmax _ 2 19 and the normal-state resistance of the 
junctions is R = 15.6 Q. After this characterization, the 
dc SQUID is operated at a given setpoint voltage Vsp us- 
ing a feedback loop that adjusts the flux via the stripline 
current HE 0. The feedback loop is used to reduce 
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FIG. 2: Frequency shift (a) and quality factor (b) plotted 
versus the normalized bias current. The voltage setpoint was 
halfway between V m in and 14iax in these measurements. 



low-frequency flux noise and flux drift and has a band- 
width of ~ 2 kHz, i.e., it does not respond to the 1 MHz 
resonator signal. 

The fundamental mode of the flexural resonator is ex- 
cited using a piezo element underneath the sample and 
the displacement of the beam is detected as follows: The 
in-plane magnetic field transduces a displacement of the 
beam u into a flux change ~ £Bu, which in turn changes 
the voltage over the dc SQUID. This voltage is amplified 
using a cryogenic high-electron mobility transistor fol- 
lowed by a room temperature amplifier and then recorded 
using a network analyzer. Figure shows the amplitude 
and phase of the measured response, from which the res- 
onance frequency /r and quality factor Q are obtained. 

To observe backaction of the dc SQUID detector on the 
resonator, the frequency response is measured for differ- 
ent bias conditions of the SQUID. Figure [2] shows that 
both fn and Q depend on the bias current /#. The 
resonance frequency saturates at fo = 1.053010 MHz for 
large positive and negative bias currents. However, when 
decreasing the /#, /r first goes up by a few hundred Hz 
around /™ ax and then it decreases rapidly with about 
-2000 Hz at the lowest stable setpoint voltage. This is 
more than 10 x the linewidth Jr/Q = 194 Hz of the reso- 
nance shown in Fig. [3J. Figure [2]o shows that the quality 
factor of the resonator changes from Qo = 5300 to less 
than 2000. Similar to the resonance frequency, first an 
increase and then a stronger decrease in Q is observed 
when lowering the bias current. Figure [3] shows that the 
frequency and damping can also be changed by adjusting 
Vsp, i.e., the flux through the SQUID loop. The shifts 




FIG. 3: Measured bias current and voltage setpoint depen- 
dence of the frequency shift (left) and damping (right). Points 
without a good lock are indicated in gray. 



are largest for low setpoints and low bias currents (dark 
regions). The regions with a lower frequency coincide 
with the regions where the damping has increased. Bias 
points with positive frequency shifts and increases in Q 
are indicated in white. Finally, by varying the driving 
power we confirm that the observed effects are not due 
to nonlinearities in the SQUID or in the resonator (iBj . 

Unlike for position detectors such as the SET where the 
backaction originates from the Coulomb force, the back- 
ward coupling between the SQUID and the beam is the 
Lorentz force Fl [l8|, |2l| . This force is due to the current 
that flows through the beam in the presence of the mag- 
netic field that couples the resonator and the SQUID. A 
displacement changes the flux, and this in turn changes 
the circulating current in the loop J [24], giving a dif- 
ferent force on the beam. In addition, resonator motion 
yields a time- varying flux through the loop, which in- 
duces an electro-motive force and thereby also generates 
currents that change the Lorentz force [26j. 

The displacement of the fundamental out-of-plane flex- 
ural mode u is given by [27]: 
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- muJoii/Qo + ttiujqU = F d (t) + F L {t). (1) 



The resonator has a mass m, (intrinsic) frequency fo = 
uo/2tt and quality factor Qq. Fd is the driving force 
and Fl = aB£(lB/2 + J) is the Lorentz force. Here, 
a = (uf) -1 Jq u(x) dx « 0.9 for the fundamental mode, 
so that also d^/du = aB£ [H, Hzl- For small ampli- 
tudes and low resonator frequencies (much smaller than 
the characteristic SQUID frequency uo c = 7ri?/™ ax /<I>o), 
the average circulating current can be expanded in the 
displacement and velocity u [25]: 



dJ 



dJ 



J(u, u) = Jq + -g^aBiu + -^-aBiu. 



d<$> 



(2) 



Inserting Eq. ((2]) into Eq. (pp) shows that the dJ/d& term 
affects the spring constant muo\ and thus whereas 
the dJ/d& term renormalizes the damping. The shifted 
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resonance frequency and quality factor are: 

Jl d2 p2 rmax 

f R = f (l-A fU ) 1/2 , with A, = — ^ ,(3) 



fo 1-A 



with An = 



mcjQ 2$ 



• (4) 



Here, j = dJ/d<£> x 2<l> // c max and ^ = dJ/d<$> x 
2cj r ^>o//^ riax are the scaled flux-to-current transfer func- 
tions [28]. The former indicates how much the circulat- 
ing current changes when the flux through the ring is 
altered, whereas the latter quantifies the effect of a time- 
dependent flux on the circulating current. These func- 
tions were first studied in the analysis of the dynamic 
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input impedance of tuned SQUID amplifiers 28, 
are intrinsic properties of the SQUID. 

Before looking in more detail at the transfer functions, 
we first focus on the coupling. The dimensionless pa- 
rameters Af [20] and Aq characterize the backaction 
strength. They contain the term aB£ squared as both 
the flux change and the Lorentz force are proportional to 
the magnetic field. This implies that the backaction re- 
mains the same when the direction of the magnetic field 
is reversed and this is what we observe experimentally. 
Af is proportional to 7™ ax , whereas the damping induced 
by the SQUID depends on R. By a careful design of the 
resonator and SQUID, the backaction strengths can be 
tuned over a wide range. Eqs. (|3j) and (|4]) show that the 
largest backaction occurs for large flux changes aB£, low 
spring constants mcj 2 ,, and large circulating currents, i.e., 
large Iq and low R. For the device studied in this Letter, 
we estimate A f = 4.1 x 10 -4 and Aq = 2.8 x 10 -4 . Fi- 
nally, note that the two coupling parameters are related 
by Aq = Af x Q uj /uj c . 

So far, the analysis did not assume anything about 
the number of junctions, nor about their microscopic de- 
tails. To obtain the transfer functions and j^, we 
model the junctions in the dc SQUID using the RCSJ 
model {24[. The transfer functions can be calculated an- 
alytically in certain limits [30|. However, to obtain their 
full bias-condition dependence, and must be cal- 
culated numerically. This is done by simulating the dy- 
namics of the SQUID in the presence of a time- varying 
flux |25j. Figure H^l shows the bias-dependence of j^. In 
the region where V = 0, the circulating current redis- 
tributes the bias current between the two junctions such 
that no voltage develops. Here the circulating current is 
of the order of 7™ ax /2, which gives ~ — 1 (blue). In 
the dissipative region (V ^ 0), the circulating current 
is suppressed. Therefore, the circulating current changes 
rapidly close to the edge of the dissipative region. The 
orange color in Fig. [4^ indicates that is large and 
positive near the critical current. The largest downward 
frequency shift is expected near a half-integer number of 
flux quanta, whereas vanishes for integer flux. With 
the value of the coupling parameter A f and the resonance 



frequency fo the frequency shift is calculated as shown in 
Fig. Ht. The maximum value = 53 gives a frequency 
shift of —12 kHz in the lower- left corner, which has to 
be compared with the experimental value of ~ — 2 kHz. 
Increasing the bias current above the critical current re- 
sults in a smaller (light yellow and light blue) that 
depends linearly on the inductive screening parameter 
Pl [24( for the experimental conditions. In this region 
the simulations predict both positive (blue) and nega- 
tive (yellow) value for j^. Positive and negative shifts 
are also observed in the experiment (Figs. [2]and[3j). The 
largest negative value found in the simulations of that re- 
gion is = — 0.65, which results in an increase in fn of 
~ 140 Hz, which is in agreement with the observed value 
of ~ 200 Hz (Fig. [2^l) . For even larger bias currents, the 
frequency shift vanishes fn ~ fo. This corresponds to 
the flat regions in Fig. [2^. 

The change in damping is determined by as indi- 
cated by Eq. (j4]). Its dependence on the bias conditions 
is shown in Fig. HJ). Well inside the experimentally in- 
accessible non-dissipative region ~ +1 and the back- 
action results in a small increase in Q. In the opposite 
limit of large bias currents = —tt (light blue). This 
value combined with the small value of Aq implies that 
the quality factor in the flat region in Fig. EJd is close to 
the intrinsic Q-factor, Qo- In this region the small addi- 
tional damping is due to the current induced by the time- 
varying flux i>/2i?, which is dissipated in the junction re- 
sistances [30| . This contribution is well-known from mag- 
netomotive readout of mechanical resonators [26[ . When 
lowering the bias current, changes sign and rises to 
about +500. This reduces the damping and might even 
lead to instability (Q < 0) if Aq is large enough. This 
decrease of damping corresponds to the bumps in Fig. 
[2)3. The largest observed quality factor Q = 5800 corre- 
sponds to = +400, which is in reasonable agreement 
with the simulations. Close to the critical current, 
goes to large negative values leading to an enhanced dis- 
sipation. Figure [4jd shows that the calculated Q-factor 
is indeed lowest near the critical current. In summary, 
our model shows that although the coupling strength is 
small, the dynamics of the dc SQUID greatly enhances 
the backaction. 

Various interesting effects can be observed when the 
backaction is strong. If the resonator and SQUID are 
strongly coupled, the resonator temperature is set by the 
effective bath temperature llJEI] of the SQUID. The in- 
creased damping cools the resonator, but the shot noise 
in the bias current leads to an increase in the force noise 
on the resonator, heating it. The question whether the 
resonator temperature is above or below the environmen- 
tal temperature should be addressed in future research. 
Furthermore, the dependence of fu and Q on the bias 
conditions allows parametric excitation of the mechani- 
cal resonator by either modulating the flux or the bias 
current. This enables squeezing of the thermomechani- 
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FIG. 4: Surface plots with iso- voltage lines at different bias 
conditions of a dc SQUID. The (logarithmic) color-scale rep- 
resents the calculated flux- to- current transfer functions (a) 
and (b). With the values for A/ and Aq, the frequency 
shift (c) and quality factor change (d) are calculated. The 
simulation is done for the experimental conditions where the 
inductive screening parameter /3l = 0.21 and the Stewart- 
McCumber parameter f3 c = 0.23 [H[25[. 

cal noise of the resonator 0, Finally, if the SQUID 
contains multiple, nearly identical mechanical resonators, 
these are coupled to each other by the backaction. This, 
in turn, can synchronize their motion and might lead to 
frequency entrainment if higher order terms in Eq. (j2j) 
become significant [32|. These examples are only a few 
intriguing possibilities of the rich physics connected to 
the backaction that we have described in this Letter. 
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Temperature and power dependence 

Figures [5^ and b show the temperature dependence of 
the resonance frequency and quality factor. These values 
are measured at large bias currents where the backaction 
is negligible. The frequency change due to temperature 
is small compared to the observed backaction (see main 
text). The intrinsic quality factor decreases significantly 
with increasing temperature. This rules out that the ob- 
served frequency shift and Q-factor change are caused 
by heating of the resonator due to Joule heating in the 
junctions: We observe an increase in quality factor with 
increasing bias current and voltage setpoint (i.e. dissi- 
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FIG. 5: Temperature dependence of the intrinsic quality fac- 
tor Qo (a) and resonator frequency fo (b). (c) Colorscale plot 
of the oscillator response at B = 100 mT at different driving 
powers Pa- This measurement is done with a voltage setpoint 
halfway between Vmin and Vmax- 



pation), but a decrease in quality factor with increasing 
temperature. An increased damping at higher temper- 
atures is seen more often in micro- or nanomechanical 
resonators [33-13^]. 

The observed frequency shift and change in damping 
do not depend on the driving power. As shown in Figure 
[5fc, the measured resonator response stays the same in all 
panels. Although the driving power is changed by three 
orders of magnitude, the only effect is that the signal-to- 
noise ratio becomes better when increasing the power. 

For the highest driving power (Pd = — 75 dBm) and 
highest Q-factor (Q ~ 5800) the amplitude of the res- 
onator motion is w ma x = 20 pm, as determined using 
the calibrated displacement responsivity [16|. The flux 
through the dc SQUID is then modulated with an am- 
plitude aB£u max = 0.02 <3>o. So even for the largest res- 
onator motion, the change in flux is much smaller than 
a flux quantum. Exactly on resonance, the piezo motion 
is Q times smaller than ii max , about 3 fm. The driving 
force Fd is then given by the resonator mass m times the 
acceleration of the piezo element, Fd = muj 2 u p « muj^Up. 
The measurements shown in the main text are done with 
a driving power of -80 dBm (Fd = 48 fN). 



Device B 

All effects that we have observed in the device that we 
discuss in the main text have also been measured in a 
second device. The resonator in this dc SQUID operates 
around 2 MHz. The maximum critical current of device 
B is J c max = 2I = 2.4 /iA at B = 115 mT and 7 c max = 
1.0 /iA at B = 130 mT. 




FIG. 6: Backaction measurements of device B. (a) Frequency 
shift and (b) quality factor plotted versus the normalized bias 
current at B — 115 mT (circles) and B — 130 mT (triangles). 
The inset depicts the voltage setpoint for these measurements 
with respect to V m in and Kiax- Bias current and setpoint 
dependence of the frequency shift (c) and damping (d) at 
B = 115 mT. The solid lines indicate the measured Vmin and 

Vmax • 



In the measurements on device B, the feedback loop 
could not maintain the SQUID voltage for low voltage 
setpoints. Also, a less sensitive, room temperature am- 
plifier was used for the resonator signal. Its lower gain re- 
sulted in an increased scatter in the data and the inability 
to explore the region with the highest backaction, i.e., at 
low bias current and low voltage setpoints. The measure- 
ments in Fig. [6] show qualitatively the same backaction 
as the data presented in the main text. 



The RCSJ model for the dc SQUID 

To calculate the backaction, the dc SQUID is modelled 
using the resistively- and capacitively-shunted junction 
(RCSJ) model. This widely-used model is discussed in 
detail in Ref. [24]. The introduction to this model pre- 
sented here is largely based on this review. A current 
Ib is sent through the SQUID and the circulating cur- 
rent J redistributes the current over the two junctions, 
which we assume to be identical. In the RCSJ model, 
the two junctions (labelled with i = 1,2) are modelled as 
a resistor (i?), capacitor (C) and an "ideal" Josephson 
junction with critical current Iq in parallel. The voltage 
over each junction is related to the time derivative of the 
phase difference Si of the superconducting wave function: 
Vi = $o Si/2n, where $ = h/2e = 2.05 x lO" 15 Tm 2 
is the flux quantum. Current conservation yields two 
second-order differential equations, governing the time- 
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dependence of the phase differences of the junctions: 

^CS 1 + -^--6 1 + losing = -Ib + J, (5) 
Ztt Ztt H Z 

p-CSi + 7^ -^2 + ^0 sin 5 2 = \l B -J. (6) 
Zir Ztt K Z 

These equations are coupled to each other by the amount 
of flux piercing the loop: 

62-6! = 2tt- $tot/*o- (7) 

The total flux <£> t ot has two contributions: the externally 
applied flux $ (which also includes the flux due to the 
resonator displacement) and the flux due to the circulat- 
ing current flowing through the inductance of the loop L, 
i.e., $tot = $ + LJ. 

The equations are scaled to simplify their numerical 
integration. This yields: 

PcSi+Si+amS! = i B /2 + 3 (8) 
Pcfa + fa + smSz = i B /2-J (9) 
27r(0 + /3 L j/2) = 6 2 -6l (10) 

The bias current and circulating current are normalized 
using the critical current: %b = Ib/Iq and J = J/Io- 
Furthermore, time is scaled using the characteristic fre- 
quency uj c = 2ttRIo/^o, fluxes using the flux quantum, 
i.e., (j) — ^/^o, and voltages using the characteristic volt- 
age Rio so that v = V/RIo = {5i-\-&2)/2. The parameter 
p L = and p c = 2kI R 2 C/<S> are the inductive 

screening parameter and the Stewart-McCumber num- 
ber respectively. The inductive screening parameter in- 
dicates how much a change in flux is screened by the 
circulating current J flowing through the self-inductance 
of the loop L, whereas Pc indicates the importance of 
inertial terms due to the junction capacitance C. 

The three equations are integrated numerically for dif- 
ferent bias conditions, i.e., different values for the bias 
current Ib and for the flux <I> through the SQUID. Fig- 
ure [7^ shows typical examples of calculated time- traces of 
the circulating current j and voltage v. Both are rapidly 
oscillating at a frequency of 0.69 co> c , which is the Joseph- 
son frequency that equals the average value of the volt- 
age v [24(. The maximum current that the dc SQUID 
can carry without generating a voltage equals the sum of 
the critical currents of the two junctions: /™ ax = 22o- 

The critical current (Iq = 1.1 /iA) and the normal-state 
resistance (R = 15.6 Q) of the junctions are estimated 
from the IV-characteristics (Fig. lb of the main text). 
Using finite-element simulations we estimate L = 175 pH 
for our device. Finally, the capacitance C = 0.6 pF is 
obtained from the position of the LC resonance in the dc 
SQUID 0. 




FIG. 7: (a) Calculated timetrace of the circulating current 
j(t) and voltage v(t) in the absence of noise, (b) Circulat- 
ing current for a small (0 mo d = 0.01) modulation of the flux 
with c^mod = 0.02. The time-averaged value of the circulating 
current j(t) has a phase shift with respect to the modulation 
<p(t) as indicated by the orange lines, (c) Absolute value of 
the Fourier transform of the timetraces of the SQUID volt- 
age and circulating current shown in (b). (d) The modula- 
tion frequency dependence of the real (dark gray) and imag- 
inary (light gray) parts of the transfer function. The lines 
are a guide to the eye. These simulations were done for the 
dc SQUID parameters from the main text (Pl = 0.21 and 
Pc 0.23) at i B = 2 and 0.25. 



Calculation of the transfer functions 

When the resonator moves, the flux through the dc 
SQUID loop is altered, which in turn changes the average 
circulating current J. In principle, J(t) could depend on 
the all the past displacements, u(t') for t' <t. However, 
the dc SQUID reacts at a frequency (~ uj c /2tt ~ 8 GHz) 
that is much faster than that of the resonator (1 MHz), so 
the circulating current J(t) is expected to depends on the 
instantaneous displacement u(t). For small amplitudes 
(u <C $0/ aB£) the response is linear and gives a contribu- 
tion J\{t) = c\u{t). Another contribution comes from the 
velocity of the resonator, ii, which causes a time- varying 
flux. This generates an electromotive force in the SQUID 
loop (Faraday's induction law), which also changes the 
circulating current by J2W = C2u{t). Combining these 
two effects gives J(t) = J{u,ii) = c\u(t) + C2ii{t). The 
values of the parameters c\ and C2 depend on the dynam- 
ics of the dc SQUID and are c\ = dJ/du and C2 = dJ/du. 
As discussed in the main text, these quantities are related 
to the intrinsic flux-to-current transfer functions, — 
<f>o/(aB£I ) x dJ/du and = u c ^ /(aB£I ) x dJ/du. 
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In principle can be obtained by calculating the 
average circulating current at different fluxes and then 
numerically differentiating this to obtain = dj/dcj). 
However, for the velocity-dependent transfer function 
this is not possible. Our method for calculating these 
transfer functions and works as follows: We cal- 
culate the steady-state response of the circulating cur- 
rent with a small modulation added to the applied flux, 
(j) —> (j) + (/>mod cos(o; mo dt). Figure [Tb shows that this 
modulates the circulating current j(t). Figure [?h shows 
the Fourier transform of the circulating current and the 



voltage. In the spectrum of both j and v a peak appears 
at the modulation frequency. The real part of the peak 
corresponds to the derivative = Re[j mo d/^mod], while 
Imbmod/^mod] = -Umodjj, and similar for v mod . The fre- 
quency dependence in Fig[7]i, shows a constant Re[j mo d] 
and a linearly increasing Im[j mo d] as indicated by the 
dotted lines. The transfer functions do not depend on 
the modulation frequency, provided that it is sufficiently 
low. This confirms that the circulating current only de- 
pends on the instantaneous displacement and velocity as 
was postulated at the beginning of this Section. 



